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Abstract 

We propose a hybrid model of portfolio credit risk where the dynam- 
ics of the underlying latent variables is governed by a one factor GARCH 
process. The distinctive feature of such processes is that the long-term 
aggregate return distributions can substantially deviate from the asymp- 
totic Gaussian limit for very long horizons. We introduce the notion of 
correlation surface as a convenient tool for comparing portfolio credit loss 
generating models and pricing synthetic CDO tranches. Analyzing alter- 
native specifications of the underlying dynamics, we conclude that the 
asymmetric models with TARCH volatility specification are the preferred 
choice for generating significant and persistent credit correlation skews. 
The characteristic dependence of the correlation skew on term to maturity 
and portfolio hazard rate in these models has a significant impact on both 
relative value analysis and risk management of CDO tranches. 

1 Introduction 

The latest advances in credit correlation modeling were in part motivated by 
the growth and sophistication of the so called correlation trading strategies, 
in particular those involving the standard tranches referencing the Dow Jones 
CDX (US) and iTraxx (Europe) broad market CDS indexes. The synthetic 
CDO market allows investors to take views on the shape of the credit loss 
distribution of the underlying collateral portfolio. The market implied portfolio 
loss distribution is now well exposed through the pricing of liquid standard 
tranches, which in turn are expressed through their implied correlations. 

The pricing of credit derivatives has been based on either structural Merton 
style models or reduced form credit migration models. See Lando [21] and 
Schonbucher [37] for a survey of these techniques. In both cases, rather ad hoc 
models of dependence are needed to explain why correlations vary over time and 
across tranches. Frequently, specific copulas are postulated to model prices at 
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a point in time. However, this approach does not easily generalize to dynamic 
situations where new information is continually being revealed and prices of 
tranches of different maturities evolve consistently with each other and with the 
underlying reference portfolio. 

This paper brings the standard time series methodology into the portfolio 
credit risk setting. Specifically, we consider a structural credit model where the 
latent asset variables evolve according to a one factor multi-variate asymmetric 
GARCH model. This model is formulated and estimated on high frequency 
data and the implications for long horizon loss distributions are derived by time 
aggregation. We demonstrate that the asymmetric GARCH specification can 
generate multivariate return distributions with both significant lower tail depen- 
dence and asymmetry that persist for very long time horizons. We then show 
that a multiname credit portfolio loss distribution derived from our model can 
produce implied correlation skews similar to the ones observed in the synthetic 
CDO market. 

The article is organized as follows. In section [2] we give a brief overview 
of the portfolio credit modeling in a general copula framework, and introduce 
the notion of correlation surface as a generalization of base correlation to a 
dynamic setting. We derive analytical formulas expressing the portfolio loss 
distribution directly in terms of the shape of the correlation surface. In section 
[3] we apply time series models to the portfolio credit risk problem. We develop 
new implications of the well known asymmetric threshold GARCH (TARCH) 
models and show their power in explaining persistent non-Gaussian features of 
long horizon market returns. We then demonstrate the implications of using 
such latent variable specification for default correlation estimates. In section [4] 
we use the correlation surface estimates to assess the ability of various static 
and dynamic credit loss-generating models to produce realistic pricing of CDO 
tranches. In section[5]we summarize the results and outline possible applications 
and extensions of our approach. 

2 Defining Portfolio Credit Risk Models 

In this section we will introduce the notion of correlation surface as a general tool 
for comparing portfolio credit risk models. We will first give a brief overview of 
the general copula framework, focusing particularly on the symmetric one factor 
latent variable assumption and the large homogeneous portfolio approximation 
which will be used throughout this paper. We will then formalize the definition 
of the correlation surface and show that, under the above approximation, it 
contains sufficient information to recover the full portfolio loss distribution and 
to price CDO tranches. 

2.1 General Copula Framework 

Consider a portfolio of M credit-risky obligors in a static setup with a fixed time 
horizon [0,T]. To simplify notations we will skip the time subscript for time 
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dependent variables. At time t = all M obligors are assumed to be in non- 
default state and at time T firm i is in default with probability pi. We assume we 
know the individual default probabilities p = \pi, ■■■,Pm] (either risk- neutral, 
e.g. inferred from default swap quotes, or actual, e.g. estimated by rating 
agencies). Let > be the random default time of obligor i and = l{ Ti <T} 
the default dummy variable which is equal to 1 if default happened before T 
and otherwise. 

The loss generated by obligor i conditional on its default is denoted as k > 0. 
The loss k is a product of the total exposure and percentage loss given default 
1 — Ri where Ri G [0, 1] is the recovery rate. We assume that all U are constant 
(see pQ for discussion on stochastic recoveries). Portfolio loss Lm at time T is 
the sum of the individual losses for the defaulted obligors 

M M 
Lm = / ] kl{n<T} = / ] kYj (1) 

i=l i=l 

The mean loss of the portfolio can be easily calculated in terms of individual 
default probabilities: 

M M 

E(L M )=Y,kE(Yi) = Y,kPi (2) 

1=1 i=l 

Risk management and pricing of derivatives contingent on the loss of the 
credit portfolio, such as CDO tranches, require knowing not only the mean but 
the whole distribution of losses with cdf Fl (x) = P (Lm < x) . Portfolio loss dis- 
tribution depends on the joint distribution of default indicators Y = [Yi, Ym]' 
and in a static setup can be conveniently modeled using the latent variables 
approach [TJ]. Particularly, to impose structure on the joint distribution of de- 
fault indicators we assume that there exists a vector of M real-valued random 
variables R = [Ri, Rm] and M dimensional vector of non-random default 
thresholds d = [d\, (1m}' such that 

Yi = 1 R, < di for i = 1, M (3) 

Denote F : R AI — > [0, 1] as a cdf of R and assume that it is a continuous 
function with marginal cdf {-Fi} i=1 . For each obligor i the default threshold di 
is calibrated to match the obligor's default probability p. t by inverting the cdf 
of its aggregate returns Ri : di — F^ 1 (pi). According to Sklar's theorem [3Tj . 
under the continuity assumption F can be uniquely decomposed into marginal 
cdfs {Fi}^ and the M-dimensional copula C : [0, 1) M => [0, 1] 

F(d)=C(F 1 (p i ),..,F M (p M )) (4) 

The most popular copula choices are the Gaussian copula model 22] , Student- 
t [23], and double-t [2U]. The choice of copula C defines the joint distribution of 
default indicators from which the portfolio loss distribution can be calculated. 
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The number of names in the portfolio can be large and therefore the calibra- 
tion of the copula parameters can be problematic. To reduce the number of 
parameters some form of symmetry is usually imposed on the distribution of 
default indicators. Gordy jTS] and Frey and McNeil [2] discuss the mathemat- 
ics behind the modeling of credit risk in homogeneous groups of obligors and 
the equivalence of the homogeneity assumption to the factor structure of default 
generating variables. 

Assumption 1 (Symmetric One Factor Model): Assume that loss given 
default li = (l-Ri)-ni and individual default probabilities pi are the same for all 
M names in the portfolio and that the latent variables admit symmetric linear 
one factor representation: 

ni — n (5a) 

Ri = R (5b) 

Pi=P (5c) 

Ri = bR m + \fl - b 2 Ei with < b < 1 (5d) 

where R m and Ei are independent zero mean, unit variance random variables. 
E[-s are identically distributed with cdf G(»). 
Within this framework, denote: 



F (di) = P (Ri < di) cdf of aggregate total returns Ri 
G(di) = P (Ei < di) cdf of aggregate idiosyncratic returns Ei 
F (d) = P (R < d) joint cdf of R 
• C (u) ee F (F- 1 (ui) , F- 1 (u M )) copula of R 

Note than the assumption of one factor structure implies that equity returns 
R are independent conditional on the market return R m and therefore F (d) 
can be computed as expectation of the product of conditional cdfs: 

F (d) = E (j[P (Ri < dilRm^J = E (j[G (di - hR m )j (6) 

Parameter b defines the pairwise correlation of latent variables with the 
market factor. The correlation of latent variables, p, which is often referred to 
as " asset correlation" (this naming reflects the interpretation of latent variables 
as asset returns in Merton-style structural default models), is constant across 
all pairs of assets in a symmetric single factor model: 

Pij =P = b 2 (7) 

In addition to the asset correlation, which reflects the co-movement of returns 
on small scale, multivariate distributions can be also characterized by measures 
that reflect joint extreme movements for a pair of assets - the tail dependence 
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coefficient Xfj and the pairwise default correlation coefficient pfj (p) . Suppose Ri 
and Rj are the stock returns for companies i and j over the [0, T] time horizon. 
The coefficient of lower tail dependence and the default correlation coefficient 
for two random variables with the same continuous marginal cdfs, F (R) , and 
the same default probabilities, p,are defined as: 

X% = lim P (Ri < d p \R 3 < d p ) = lim (8) 

p—y+O p-v+0 p 

d t \ /, , \ C{p,p)-p 2 

Pij (P) = corr{l {Ri < dp} ,l {Rj < dp} ) = _ (9) 

where p is the probability of crossing the threshold (also interpreted as the de- 
fault probability), and is related to the latter via the relationship d p = F^ 1 (p) . 
Both these measures depend only on the bivariate copula of the two random 

variables and are asymptotically equal: lim pf- (p) = Af -. Generally speaking, 

p— H-0 

the measures of small-scale and extreme co-movement of assets are independent 
of each other. For example, it is quite possible to have p^ — and Xfj ^ 
and vice versa for a non-Gaussian multi-variate distribution. Embrecht et al 
provide very detailed introduction to the properties of those dependence 
measures. 

To simplify the calculations even more, the large homogenous portfolio (LHP) 
approximation is often used. Suppose that we increase the number of names 
in the portfolio while keeping the total exposure size of the portfolio con- 
stant so that rii — N/M. Conditional on R m the loss of the portfolio con- 
tains the mean of independent identically distributed random variables, Lm — 
^1 — Rj Njj X)f=i l{.Ri<d}7 which a.s. converges to its conditional expectation 
as M increases to infinity. We use L without subscript to denote the portfolio 
loss under LHP assumption. 

Proposition 1 (LHP Loss) Under Assumption 1 



L = lim 

M->oo 



M 

" {Ri<d} 



i=l 

= ^1 — Rj NG ( ^ a.s. for any R m £ supp(G) 

Proof, see proposition 4.5 in [14] ■ 



1-R )NP(Ri<d\R m ) (10) 



Based on ( 10 ) cdf of L can be expressed in terms of the cdf of R m 

P{L<l)=P{R m >d 1 {l)) (11) 



We use the following notation for the Gaussian distribution: 
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Notation 2 Let $ (.) and (/>(.) with one agrument denote the cdf and pdf, cor- 
respondingly, of a standard normal random variable. Let $ (., .; p) and </>(., .; p) 
with three arguments denote cdf and pdf of two standard normal random vari- 
ables with linear correlation p, and $j (., .; p) denote the partial derivative of 
$ (., .; p) with respect to the i'th agrument, e.g. $3 (., .; p) = (., .; p) . 

For the Gaussian copula we have the familiar formula for LHP loss first de- 
rived by Vasicek [33J , where we have substituted the asset correlation parameter 
p in place of the factor loading b using the relation Q: 

L G = (l - R) iV$ ( ^ l{P ^ Rm ) (13) 
P(i<0 = l-«(df (0) (14) 

tfW-^-^W^-l (15) 

Vp Vp \ i-r)n 



Vasicek [33] and Schonbucher and Shubert [28] show that LHP approxima- 
tion is quite accurate for upper tail of the loss distribution even for mid-sized 
portfolios of about 100 names. We will use a symmetric one factor LHP ap- 
proximation in this paper for analytical tractability. 

2.2 Prom Loss Distribution to Correlation Surface 

It is intuitively clear that the choice of the dependence structure affects the 
degree of uncertainty about the portfolio loss. Indeed, if all issuers in the port- 
folio are completely independent of each other and any common driving factor, 
then the law of large numbers assures that the portfolio loss under the LHP 
approximation is a well determined number with little uncertainty about it. On 
the other extreme, if the issuers are highly dependent such that they all default 
or survive at the same time, then the portfolio loss has a binary outcome - it is 
either zero or equal to the maximum loss. 

It is easy to quantify this intuitive result. While the mean of the loss distri- 
bution is not affected by the choice of copula, one can show that the second and 
higher moments of the loss distribution depend on the copula characteristics. In 
particular, the variance of the loss can be expressed in terms of bivariate default 



correlation coefficient p d (p) defined in section 2.1 Under assumption of equal 



default probabilities for all obligors, it is given by: 

Var{L) = (1 - R) 2 N 2 p{\ - p)p d (p) (16) 

Thus, in line with the intuition, the uncertainty of the default loss distribution 
is directly proportional to the default correlation coefficient. 

To characterize the shape of the entire distribution of L, we must look at the 
particular slices of portfolio loss. Let (Kd,K u ] denote a tranche with attach- 
ment point Kd and detachment point K u expressed as fractions of the reference 
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portfolio notional so that < Kd < K u < 1. The notional of the tranche is 
defined as N(x d ,K u ] = N (K u — Kd) where N is the notional of the portfolio. 
The loss L(x d ,K u \ of the tranche is the fraction of L that falls between Kd and 
K u . For simplicity, assume that total notional N is normalized to 1. 



L(K d ,K u ] = f(K d ,K u ] (L) (17) 

f ( K d ,K u ] (x) = (x- K d y - (x - K u )+ (18) 

Tranches with zero attachment point, (0, K u ] , and unit detachment point, 
(Kd, 1] , are called equity and senior tranches, respectively. Loss of any tranche 
can be decomposed into losses of either two equity or two senior tranches 
tranches L( KdiKu ] = L( ,k„] - L (o,K d ] = L (K d ,i] - l (k u ,i] This is similar to 
representing a spread option as a long/short position using either calls or puts. 

The expected loss of the equity tranche i(o,if] depends on the portfolio loss 
distribution and under the LHP approximation can be computed using only the 
distribution of the market factor: 



= {1-R)E 



G 



d - bR m 



{R m >ch(K)} 



(19) 

KP [R m < d x (K)) 



The expectation in ( 19 ) can be computed by Monte Carlo simulation or numer- 



ical integration if we know the cdf of residuals G and the distribution of R m 
(see appendix [C]) . For the Gaussian copula, the integral can be calculated in a 
closed form: 



E G L {QtK] = (l - k) $ 0) , -di; -y/p) + {di) (20) 

dr - ^ (P) - (^J) (21) 

Because of its analytical tractability, it is convenient to use the Gaussian 
copula as a benchmark model when comparing different choices of dependence 
structure. By finding the asset correlation level p that replicates the results 
of more complex portfolio loss generating models in the context of a Gaussian 
copula framework, we can translate the salient features of such models into 
mutually comparable units. More specifically, we define the correlation surface 
as follows: 

Definition 3 Suppose the loss distribution of a large homogeneous portfolio is 
generated by a model ic,p, i?| with copula C, identical individual default proba- 
bilities p and recovery rates R. Let L(n t K\ G 



P f (0> K] f(0,K] ((l-R)p) 
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be the expected loss of the equity tranche (0, K] . We define the correlation sur- 
face p(K,p, R) of the model |c,p, J?| as the correlation parameter of the Gaus- 
sian copula that produces the same expected loss EL^x] for the tranche (0, K] 
for the given horizon T and given single-issuer cumulative default probability p: 

p(K,p, R) solves E G L (0 , K] (p) = EL {Q , K] for all K e [0, 1] (22) 
where EL^ x] * s expected loss of the tranche 

(0,K] generated by model ^C,p,R^ 



where E Liq^k] *s defined in (20) 



The correlation surface as denned above is closely related but not identi- 
cal to the notion of the base correlation used by many practitioners [2 5) . The 
difference is that the base correlation is defined using the prices of the equity 
tranches, which in turn depend on interest rates, term structure of losses, etc. 
By contrast, the correlation surface is defined without a reference to any mar- 
ket price. It characterizes the portfolio loss generating model, rather than the 
supply /demand forces in the market. 

One could, of course, take another logical step, and instead of deriving the 
correlation surface from the parameters of the dynamic loss generating model, 
go in the opposite direction - find such parameters of the loss generating model 
that result in the closest match to the market prices of CDO tranches. It would 
be natural to call this solution " implied parameters" , and the corresponding 
function p(K,p,R) "implied correlation surface". The latter would, in fact, 
coincide with the conventionally defined base correlation. 

To ensure that the correlation surface is well defined we need to prove that 



( 22 ) has a unique solution. Let us first prove the following: 



Proposition 4 For the Gaussian copula, the expected loss of an equity tranche, 
E g L(q K ^, is a monotonically decreasing function of p and attains its maximum 
when p is equal to and minimum when p is 1. 

Proof. Using Notation 2, Eq. ( 20) and the properties of Gaussian distributioi^ 



1 The following properties of two dimensional Gaussian cdf are used in the calculation 

d fx — on 

*2 0, y; p) = — $ 0, y; p) = 4> (y) *' 



8y 

d 

*3 (x,y;p) = —$(x,y;p) = <f>(x,y;p) 
op 

First formula is derived by taking the derivative and re-arranging the terms. The proof of 
the second can be found in Vasicek(l34|'l 
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we derive 

d 

EpL( ,K] = ^E G L (0 ^ K] (23) 



(i -A) 



^=4> 3 (i- 1 (p) , -d i; - VP) + 4>2 (p) , -d i; - VP) ^di 



+ A"0(di)^di 

= -^|0($- 1 (rf,-di;-Vp)<O (24) 

for any p £ (0,1). ■ 

Therefore, there is a one-to-one mapping between loss distribution and cor- 
relation surface, and our transformation does not lead to any loss of information. 
The next proposition shows how to calculate the loss cdf using the correlation 
surface and its slope along the Zf-dimention. 

Proposition 5 Suppose p(K,p, R) is the correlation surface for model |c,p, iij 
and the probability distribution function of the portfolio loss is a continuous 
function then the loss cdf can be computed from the correlation surface: 

P(L<K) = P G (L<K) + p K E G L^ K] (25) 

where p = p(K,p,R) is the level of the correlation surface, pk = -^p(K,p,R) 
is the correlation surface slope and 

P G (L < K) = 1 - $ (di) (26) 

E G L [Q . K] = {l-R) ^ (Z- 1 (p) , -dv-y/p) (27) 

dl = (p) - ( JL) (28) 

VP VP \i-rJ 

Proof, first note that the derivative with respect to K of the expected tranche 's 
loss under true copula C is related to the cdf of the loss 

^ EL (o,k ] = ^E(L-(L-K) + ) (29) 
= --^E (L - K) + = E1 {L _ K > 0} = 1 - P (L < K) 
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therefore 



P(L<K) = 1-^EL^ K] (30) 

= 1 - $ (di) + ^— ^ l*- 1 (p) , -d i; P k 

= P G (L < K) + (l - i?) ^0 (p) , -di; -VP) PA' 



where partial derivative with respect to K is computed as 
d 

dK 



e k l (o,K] = ^ eGl (o.,k] (31) 



1 - R) $ 2 (p) , -d i; -y/p) ~d x + K4> (d x ) ~dx + $ (di) 



$(di 



Another important point is that the correlation surface depends implicitly 
on the term to maturity via the cumulative default probability p. However, this 
is not the only dependence - potentially the dependence structure characterized 
by the copula C also exhibits some time dependence when viewed within the 
context of the Gaussian copula. This statement needs a clarification - the copula 
C itself is defined in a manner that encompasses all time horizons and therefore 
cannot depend on any particular horizon. However, when we translate the 
tranche loss generated with this dependence structure into the simpler Gaussian 
model the transformation that is required may depend on the horizon T. In 
other words, the shape of the correlation surface p(K,p, R) may depend on the 
horizon. We study this dependence in detail in section [4| 

These results allow us to further develop the often mentioned analogy be- 
tween the role that default correlation plays in tranche pricing on one hand and 
the role that implied volatility plays in equity derivatives pricing, on the other. 



The starting point of this analogy is the result (16), which shows that default 
correlation is directly related to the uncertainty of portfolio loss. If we recall 
that the CDO equity tranche is essentially a call option on portfolio survival, it 
becomes clear that the price of the equity tranche should be positively related to 
the default (asset) correlation, just as the price of any equity option is positively 
related to the implied volatility of its underlying stock. 

The correlation surface takes this analogy one step further. Just as the 
implied volatility surface is sufficient to derive the risk-neutral distribution of 
the underlying stock (see [8], [9], [30], [31]) and price any European option, we 
have shown that the correlation surface is sufficient to derive the full portfolio 
loss distribution and price any CDO tranche. 
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3 Time Series Approach to Tail Risk 



3.1 Motivation for the time series approach 

The analogy between the correlation surface and implied volatility surface in- 
troduced in the previous section leads to further insight about the origins of 
the large credit portfolio loss risks. Both equity and credit derivatives pricing 
exhibits substantial deviations from the simplest Gaussian models of the under- 
lying assets. In particular, the equity index options implied volatility exhibits a 
steep downward skew of implied volatility, and the CDO tranches exhibit a steep 
upward skew of implied (base) correlations. To further underline the similarities 
between these skews, note that a senior CDO tranche with an attachment point 
that is higher than the expected loss on the underlying portfolio can be thought 
as an out-of-money put option on portfolio losses. When looked from this angle, 
the equity implied volatility skew and credit correlation skew are tilted the same 
way - towards the farther out-of-money options. 

Recall that the empirical distribution of returns does indeed exhibit signif- 
icant downside tails, and that a large part of the implied volatility skew can 
be explained by the properties of the empirical distribution [7 . Given the 
above mentioned analogies between the synthetic CDO tranches and equity in- 
dex options, it is quite natural to look for a similar explanation of the implied 
correlation skew. 

The standard Gaussian copula framework implicitly relies on the Merton- 
style structural model for definition of default correlations. Therefore, if we are 
to give an empirical explanation to the observed base correlation skew we must 
start by giving an empirical meaning to the variables in this model. Our working 
hypothesis in this paper will be that the meaning of the "market factor" in the 
factor copula framework is the same as the market factor used in the equity 
return modeling. As such, it is often possible to use an observable broad market 
index such as S&P 500 as a proxy for the economic market factor, with an added 
convenience that there exists a long historical dataset for its returns and a rich 
set of equity options data from which one can glean independent information 
about their implied return distribution. 

This hypothesis is not uncommon in portfolio credit risk modeling - for 
example, the authors of emphasized the importance of using a fat-tailed 
distribution of asset returns in the copula framework in part by citing the em- 
pirical evidence from equity markets. However, most researchers have focused 
on the single-period return distribution characteristics. 

In contrast, we focus in this paper on the long-run cumulative returns, and 
prove that their distribution is quite distinct from that of the short-term (single- 
period) returns. As we will show in the rest of this paper, it is the time aggre- 
gation properties and the compounding of the asymmetric volatility responses 
that make it possible to explain the credit correlation skew for 5- or even 10-year 
horizons. Moreover, this dynamic explanation of the skew allows one to make 
rather specific predictions for the dependence of this skew on both the term to 
maturity and on the hazard rates and other model parameters. 
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We first specify time series properties of stock returns for high frequency 
time intervals (daily or weekly) and then derive the distribution of stock prices 
over longer horizons measured in months or even years. Assuming geometric 
Brownian motion of stock prices on short intervals leads to the same log-normal 
shape for the distribution of stock prices for all future horizons. Models with 
more realistic dynamics can lead to richer distributions of the time aggregated 
returns even if the high frequency shocks are Gaussian. 

In the GARCH family of models [12] , [3] , [6] , there have been many investi- 
gations of the difference between the conditional and the unconditional distri- 
butions. These models reveal an important explanation for the excess kurtosis 
in financial returns but generally show no reason to expect skewness in returns. 
Dicbold [3J investigates the implications of time aggregation of GARCH models 
concluding that eventually they lose their excess kurtosis as the central limit 
theorem leaves the average distribution Gaussian. 

Asymmetric volatility models were introduced by Nelson |26j and further 
investigated by Glosten, Jaganathan and Runkle [TS] and by Zakoian [33]. The 
model to be discussed in this paper is essentially the GJR model and will be 
called the Threshold-ARCH or TARCH. These models all show that negative 
returns forecast higher volatility than positive returns of the same magnitude. 
This observation is very widespread and is sometimes called the leverage effect 
following Black 2 . However, it is most likely due to risk aversion as in Campbell 
and Hentschell [1] and in this paper we will refer to it simply as asymmetric 
volatility. 

The first indication that skewness could arise from time aggregation was pre- 
sented in Engle [13] . He showed that with asymmetric volatility, the skewness of 
time aggregated returns could be more negative than the skewness of the indi- 
vidual innovations. In the next section we show analytically how this negative 
skewness depends upon time aggregation. We then examine the multivariate 
distribution when the single common factor has TARCH dynamics. 

3.2 Univariate model: TARCH(1,1) 

Let Tt be the log-return of a particular stock or an index such as SP500 from 
time t — 1 to time t . Ft denotes the information set containing realized values 
of all the relevant variables up to time t. We will use the expectation sign with 
subscript t to denote the expectation conditional on time t information set: 
E t {.) = E (.\F t) ■ The time step that we use in the empirical part is 1 day or 1 
week. Predictability of stock returns is negligible over such time horizons and 
therefore we assume that the conditional mean is constant and equal to zero: 

m t = E t _ 1 (r t )=0 (32) 

The conditional volatility of = E t -i(rf) of r t in TARCH(1,1) has the au- 
toregressive functional form similar to the standard GARCH(1,1) but with an 
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additional asymmetric term [T5], [33] : 

n = otBt (33) 
a 2 = u + cirl_ x + otdr^lfa^oy + 

where {e t } are iid , have zero mean, variance normalized to 1, finite skewness 
s £ and finite kurtosis k e . We also assume that u> > and a, a d , /3 are non- 
negative so that the conditional variance of is guaranteed to be positive. 

The persistence of volatility in the model is governed by the parameter £ : 

C = E (p + as 2 + a d e 2 l {£t < 0} ) = (3 + a + a d v d e (34) 

where v d = E (e| l{ £t <o}) is the "right truncated variance" of e t ■ If C € [0, 1) 
then conditional variance mean-reverts to its unconditional level a 2 = E (<J 2 ) — 
jz^- The following parameter £ will also be useful in describing the higher 
moments of TARCH(1,1) returns and volatilities: 



£ = E {(3 + ae\ 



where k d 



2 , \2 

ad£ t l{ £t <o}J 



/3 2 



Y 2 K 



+ a 2 d k d e + 2a/3 + 2a d /3v* + 2aa d k d E 

(35) 

E {ef l{ £t <o}) is the "right truncated kurtosis" of e t . We can 



rewrite (33 1 in terms of the increments of the conditional volatility Acrf +1 



't+i 



ofand the volatility shocks rjt 



r t = ottt (36) 
A^ +1 = (1 - C) {a 2 - a 2 ) + al Vt 

T)t = a {e 2 t -l)+a d {e 2 t l {et < 0} - vf) 

The speed of mean reversion in volatility is 1 — C an d is small when £ is close 
to one which is usually true for daily and weekly equity returns - hence the 
persistence of the volatility. The TARCH(1,1) volatility shocks rj t are iid, with 
zero mean and constant variance var(rj t ) = var(ae 2 + a d e 2 l{ £t <o}) = £ — Q 2 ■ 
The correlation of conditional volatility with the return in the previous period 
depends on the covariance of return and volatility innovations: 

corr t _ 1 (r t , a 2 t+1 ) = corr t _ x (e t , r} t ) = QS g + adS e ( 37 ) 

V? - C- 

where s d = E (e| l{ £t <o}) < is the "right truncated" skewness of £t.The 
negative correlation of return and volatility shocks, often cited as the "leverage 
effect'f] is the main source of the asymmetry in the return distribution. We 



can see from formula (37) that negative return- volatility correlation can be 



2 Though we note here that the magnitudfc of this "leverage effect" in return time series 
for stocks of most investment grade issuers far exceeds the amount that would be reasonable 
based purely on their capital structure leverage. 



13 



achieved either through negative skewness of return innovations s e < 0, through 
asymmetry in volatility process a d > or combination of the two. We call these 
static and dynamic asymmetry, respectively. 

In this paper we are interested in the effects of the volatility dynamics on the 
distribution of long horizon returns which in the log representation is a sum of 

t+T 

short term log returns R^t+T = lnS , f+ x — ln^t = r u . While a closed form 

u=t+l 

solution for the probability density function of TARCH(1,1) aggregated returns 
is not available, we can still derive some analytical results for its conditional 
and unconditional moments: volatility, skewness and kurtosis. Since TARCH 
is linear autoregressive volatility process, the conditional variance Vt,t+T of the 
log return Rt,t+T is linear in <t 2 +1 : 



V t , t+T = E t Rl t+T =Et [ J2 o* = T (a 2 + (cx 2 +1 - a 2 ) 1 

(38) 

The new result of this paper for the TARCH(1,1) model is the representation 
of the skewness term structure derived in the following proposition. 

Proposition 6 Suppose < C < 1 and the return innovations have finite skew- 
ness, s £ , and finite "truncated" third moment, sf. Then the conditional third 
moment of Rt,t+T has the following representation for TARCH(1,1) 

E t Rl t+T = s E Y,E t (af + J + 3 (as e + a d s d e ) ]T ^ E t (of + u ) (39) 

u=l u=l ' 

In addition, if Ea\ is finite, then unconditional skewness of Rt,t+T is given by 



c _ H,t+T 
JT = 



ER **^ r i , Q i f , d , t(i-o-i + c t 



E(Rl t+T ) 3/ 



T l/2" e ' 3 y3/2 (^ e +a d S e ) (1-C) 2 



E 



a 

(40) 

Proof. See appendix pf] for the details. ■ 

The conditional third moment is a function of the conditional term structure 
of of, term horizon T and volatility parameters. The conditional skewness can 
be computed using second and third conditional moments derived above. The 
asymmetry in the return distribution arises from two sources - skewness of re- 
turn innovations and asymmetry of the volatility process. Note that the second 
term in the formulas for conditional and unconditional skewness is directly re- 
lated to the correlation of return and volatility innovations. If return-volatility 
correlation is zero (as e + a d s^ = 0) then St — ifT/2 s sE f^M . If return inno- 
vations are symmetric then asymmetric volatility drives the asymmetry in the 
return distribution. In figure [T] we show conditional and unconditional skew- 
ness term structures. For realistic parameters corresponding approximately to 
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Conditional Skewness Term Structure 

, 




-1 1 1 1 1 1 1 

100 200 300 400 500 600 
T(w eeks) 

Figure 1: Term structure of conditional skewness of time aggregated return 
i?t+i i t+7'.TARCH(l,l) has persistence coefficient £ = 0.98 and the following 
parametrization: a 2 = 1, a = 0.01, a d = 0.10, /3 = 0.92, s t ~ iV(0,l). We 
plotted unconditional skewness term structure and conditional for three different ini- 
tial volatilities: c 2 /2, cr 2 and 2a 2 . The term structure of Et<rf, u w&s computed from 
10,000 independent simulations. 



parameters of the TARCH(1,1) estimated for weekly SP500 log returns, both 
conditional and unconditional skewness is negative. It decreases in the medium 
term, attains the minimum at approximately the 2 year point and then decays 
to zero as T increases. The skewness term structure conditional on the high/low 
current volatility is above/below the unconditional skewness. 

To provide some empirical context to the theoretical discussion above, let 
us consider the time series of SP500 returns. The results of the estimation of 
various TARCH(1,1) specification are shown in the appendix |B| Figure [2] shows 
the estimate of skewness for overlapping returns of different aggregation horizons 
measured in days. The full sample shows high negative skewness for one day 
return because of the 1987 crash. On the post 1990 sample negative skewness 
rises with aggregation horizon up to 1 year and then slowly decays toward zero. 
Both samples show significant skewness for horizons of several years. We should 
note that confidence bounds around skewness curves are quite wide due to the 
persistence and high volatility of the squared returns and serial correlation of 
the overlapping observations. Nevertheless, both the shape and the level of the 
empirical skewness in figure [2] and the theoretical estimate shown in figure [T] are 
quite similar. 
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Figure 2: Term structure of skewness for SP500 time aggregated log returns estimated 
with overlapping samples moments for full and post-1990 data. 



3.3 Multivariate model with TARCH(1,1) factor volatility 
dynamics 

Let us now turn to a multi-variate model of equity returns for M companies, 
with a simple dynamic factor structure decomposing the returns into a common 
(market) and idiosyncratic components. To concentrate on the time dimension 
of the model we assume a homogeneity of cross-sectional return properties, 
namely that factor loadings and volatilities of idiosyncratic terms are constant 
and identical for all stocks. Thus, our homogeneous one factor ARCH model 
has the following form. 

n,t = br m . t + ae i<t 

where 

• b > is the constant market factor loading and it is the same for all stocks 

• r m i is the market factor return with zero conditional mean Et—i(r mj t) = 0, 
conditional volatility of n t = £' t _i(r^ 1 1 ) that has TARCH(1,1) parametriza- 
tion (33) 

• aEi t t are the idiosyncratic return components with constant volatilities a 2 
and zero conditional means E t -\(<J£i y t) = 

• {£i,t,Sm,t} are unit variance iid shocks for each t and all i 



Because of the simple linear factor structure and constant market loadings 

T 



time aggregated equity returns Ri t x = / Jj,u also have a one factor represen- 
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tatior0 



where R m ,T 

3.4 Tail risk and default correlation estimates 

Assuming that the latent variables in the copula framework follow the one factor 
TARCH (1,1) dynamics, we can calculate the coefficient of lower tail dependence 
Xfj and the default correlation coefficient pf j (p), which were defined in equa- 
tions ([8| and ([9]), respectively. Default correlation coefficient for Gaussian and 
Student-t can be computed in closed form, while the factor GARCH/TARCH 
models require Monte Carlo simulation, which is desribed in appendix C. Figure 
shows the numerical estimates of the default correlation pf 2 as a function of p 
for 4 different models - TARCH, GARCH and Student-t and Gaussian copulae. 

The linear correlation of latent returns is set to 0.3 for all 3 models. TARCH 
and GARCH are calibrated to have volatility dynamics parameters correspond- 
ing approximately to the weekly SP500 returns and the time aggregation horizon 
is set to 5 years. The degrees of freedom parameter for the Student-t idiosyn- 
cracies and T-Copula is set to be equal to 12. 

The GARCH distribution is symmetric and has smaller tail dependence for 
both upper and lower tails. We can see on the graph that it also has lowest 
default correlation for all default probabilities in the range of [0.01,0.2]. The 
Student-t copula is also symmetric but has fatter joint tails compared to the 
GARCH. Its default correlation is above GARCH for all p and converges to a 
positive number (the tail dependence coefficient) as p decreases to zero. 

We can see that TARCH has higher default correlation than other 2 models 
and is upward sloping for very low quantiles. The upturn for the extreme 
tails is a consequence of the left tail shape of the common factor. The default 
correlation for very low default probabilities should be close to 1 since the left 
tail of the factor is fatter than the left tail of the idiosyncratic shocks. As 
we showed in the previous sections, both kurtosis and skewness of the market 
factor declines faster for GARCH than TARCH given the same level of volatility 
persistence. 

We also show the 95% confidence bounds for default correlation, which are 
calculated using 1000 independent repetitions of 10,000 Monte Carlo simulations 
of the common factor. TARCH, T-Copula and GARCH models have virtually 
non-overlapping confidence bounds for pf • (p) for sufficiently low values of the 
default probability [0.01,0.1]. Since the correlation skew for sufficiently large 
values of the detachment point K is also related to the limit of low default 
probabilities, we can therefore rely on this result in the next section to assert 
that the differences between the correlation surfaces generated by these three 
models are statistically significant. 

3 To simplify the notations we assume that the initial time t = and use only subscipt for 
the time aggregation horizon T. 
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Ri,T = bR„ hT + Ei, T (41) 

T T 

V^?*ro,u an d ^i,T = cr^ ^£j,u are independent conditional on F o- 



u=l 




Gaussian 
Student-t 
^ — GARCH 
— TARCH 



0.02 0.04 0.06 0.08 



0.12 0.14 0.16 0.18 



Figure 3: Default correlation as a function of p for Gaussian cop- 
ula, T-Copula, GARCH and TARCH models. The linear correlation 
parameter is 0.3 for all four models. TARCH(GARCH) parameters 

are a = 0.01(0.06), a d = 0.1(0), (3 = 0.92(0.92). The degrees of freedom 
parameter for T-Copula and TARCH Student-t idiosyncracies is 12. 



4 Comparing Credit Portfolio Loss Generating 
Models 

Our goal in this paper is to provide a general framework for judging the versa- 
tility of various portfolio credit risk models. All such models, whether defined 
via dynamic multivariate returns model like in this paper or in various versions 
of the static copula framework (pQ, [T3], [IB], [25], [53] and [IE])) can be char- 
acterized by the full term structure of loss distributions. Thus, without loss of 
generality, we can refer to all models of credit risk as loss generating models, 
with an implicit assumption that any two models that produce identical loss 
distributions for all terms to maturity are considered to be equivalent. 

The correlation surface, introduced in section |2.2[ conveniently transforms 
specific choice of a loss generating model into a two dimensional surface p(K, T) 
of the Gaussian copula correlation parameter, with the main dimensions being 
the loss threshold (detachment level) K and the term to maturity T. All other 
inputs such as the recovery rate R, the term structure of (static) hazard rates 
h, the level of linear asset correlation the Student-t degrees of freedom i/, 
various GARCH model coefficients, etc. - are considered as model parameters 
upon which the two-dimensional correlation surface itself depends. 

Note that in the previous sections we have expressed the correlation surface 
as a function of detachment level and the underlying portfolio's cumulative ex- 
pected default probability p rather than the term to maturity T . Given our 
assumption of the static term structure of the hazard rates h these two formu- 
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Figure 4: Correlation surface slices corresponding to 1-year (left chart) and 5-year 
(right chart) horizons with default probabilities 0.02 and 1 — (1 — 0.02) 5 = 0.0961 for 
Gaussian, Student-t Copula with 12 degrees of freedom, and Double-t Copulas with 
(v m =12, Vj =100) and (v m =12, Vj =12) degrees of freedom. Linear correlation is 
0.3 for all copulas. 



lations are equivalent. In this section we prefer to emphasize the dependence 
on maturity horizon in order to facilitate the comparison with base correla- 
tion models and also to analyze the dependence on the level of hazard rates 
separately from the term to maturity dimension. 

Of course, the correlation surface of a static Gaussian copula model [35] is a 
flat surface with constant correlation across both detachment level K and term 
to maturity T. Any deviation from a flat surface is therefore an indication of a 
non-trivial loss generating model, and we can judge which features of the model 
are the important ones by examining how strong a deviation from flatness they 
lead to. 



4.1 Models with static dependence structure 

Let us begin with the analysis of one of the popular static loss generation models. 
On figure [4] we show the correlation surface computed for the Student-t copula 
with linear correlation p = 0.3 and v = 12 degrees of freedom. Student-t 
copula is in the same elliptic family as the Gaussian copula but has non-zero 
tail dependence governed by the degrees of freedom parameter. As a model of 
single-period asset returns the Student-t distribution has been shown to provide 
a significantly better fit to observations than the standard normal [23 . 

However, from the figure [4] we can see that the static Student-t copula does 
not generate a notable skew in the direction of detachment level K, and in fact 
generates a mild downward sloping skew for very short terms, which is contrary 
to what is observed in the market. The main reason for this is the rigid structure 
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of this model, with the tails of the idiosyncratic returns tied closely to the tails of 
the market factor. This follows from the representation of the Student-t copula 
as a mixture model [29] , Instead of producing a varying degree of correlation 
depending on the default threshold, the Student-t copula model simply produces 
a higher overall level of correlation. 

On the other hand, the more flexible double-t copula model [30] produces 
a steep upward sloping skew, as can be seen from figure [4] The main feature 
of the double-t copula that is responsible for the skew is the cleaner separation 
between the common factor and idiosyncratic returns - there is no longer a 
single mixing variable which ties the two sources of risk together. As a result, 
the idiosyncratic returns get efficiently diversified in the LHP framework and 
their contribution becomes progressively smaller for farther downside returns. 
Since the higher value of the detachment level K corresponds to farther downside 
tails, the greater dominance of the market factor translates into higher effective 
correlation for higher K, i.e. upward sloping correlation skew. 

Furthermore, by making the fully independent idiosyncratic returns more fat 
tailed one achieves a steeper skew - compare the two examples of the double- 
t copula, with the degrees of freedom of the idiosyncratic returns set to 100 
(i.e. nearly Gaussian case) and to 12 (i.e. strongly fat-tailed case), respectively. 
Indeed, for the same detachment level K the idiosyncratic returns with lower 
degrees of freedom (stronger fat tails) are less dominated by the market factor, 
resulting in relatively lower effective correlation. Since the difference between 
two cases diminishes as K grows, this translates into steeper correlation skew 
for fatter-tailed idiosyncratic returns. 

Finally, we observe that the slope of the correlation skew gets flatter as the 
time horizon grows. Within the context of double-t copula this is simply because 
the same detachment level K corresponds to less extreme tails when the term 
to maturity is greater. Following the same logic as above, this means less steep 
correlation skew. 

All of these features will have their close counterparts in the dynamic models 
which we will consider next. 

4.2 Multi-period (dynamic) loss generating models 

Let us now turn to loss generating models based on latent variables with multi- 
period dynamics. We have concluded in the previous section that a clean sep- 
aration of the market factor and the idiosyncratic returns appears to be a pre- 
requisite for producing an upward sloping correlation skew. Fortunately, the 
dynamic multi-variate models which we considered in section [3] all have this 
property, both for single-period and for aggregated returns. 

On figure [5] we show the correlation surface computed for a loss generating 
model based on GARCH dynamics with Gaussian residuals, with a linear cor- 
relation set to p — 0.3, and GARCH model parameters taken from the weekly 
SP500 estimates in appendix |B| The loss distributions and the correlation sur- 
faces are calculated using a Monte Carlo simulation with 100,000 trials. 
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Figure 5: Correlation surface for GARCH model (a=0.045, /3=0.948) with Gaussian 
shocks and the slices of the surface for 1, 3, 5 and 7 year maturities. 



As we can see, this model does exhibit a visible deviation from the flat cor- 
relation surface for short maturities. However, as we already noted in section [3j 
the distribution of aggregate returns for the symmetric GARCH model quickly 
converges to normal. Therefore, it is not surprising to see that the correlation 
surface also flattens out fairly quickly and becomes virtually indistinguishable 
from a Gaussian copula for maturities beyond 5 years. Thus, we conclude that 
the symmetric GARCH model with Gaussian residuals is inadequate for descrip- 
tion of liquid tranche markets where one routinely observes steep correlation 
skews at maturities as long as 7 and 10 years. 

Based on the empirical results of Section 3.3 we know that a GARCH model 
with Student-t residuals provides a better fit to historical time series of equity 
returns. A natural question is whether allowing for such volatility dynamics can 
lead to a persistent correlation skew commensurate with the levels observed in 
synthetic CDO markets. 

The results of section [373] suggest that the additional kurtosis of the single- 
period returns represented by the Student-t residuals does not matter very much 
for aggregate return distributions at sufficiently long time horizons. Indeed, 
figure [6] shows that the GARCH model with Student-t residuals exhibits a cor- 
relation skew that is quite a bit steeper at the short maturities, yet is almost 
as flat and featureless at the long maturities as its non-fat-tailed counterpart - 
there is a small amount of skew at 10 years, but it is too small compared to the 
steepness observed in the liquid tranche markets. Thus, we conclude that one 
has to focus on the dynamic features of the market factor process in order to 
achieve the desired correlation skew effect. 

Our next candidates are the TARCH models with either Gaussian or Student- 
t return innovations. We have seen in section [3] that the asymmetric volatility 
dynamics of these models leads to a much more persistent skewness and kur- 
tosis of aggregated equity returns that actually grow rather than decay at very 
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short horizons, and survive for as long as 10 years for the range of parameters 
corresponding to the post-1990 sample of SP500 weekly log-returns. Hence, our 
hypothesis is that a latent variable model with TARCH market dynamics might 
be capable of producing a non-trivial credit correlation skew for maturities of 
up to 10 years. 

Figures [7] and [8] show the correlation surfaces for the TARCH-based loss 
generating models. The most immediate observation is that both versions of the 
model produce a rather persistent correlation skew. Although the correlation 
surface flattens out with growing term to maturity, the steepness of the skew 
is still quite significant even at 10 years. Just as in the case of the symmetric 
GARCH model, the fat-tailed residuals lead only to marginal steepening of the 
correlation surface compared to the case with Gaussian residuals. 

Contrast these properties of the dynamic GARCH-based models with the 
features of the static double-t copula. Upon a closer inspection of figures [4] and 
[3 we can see that the TARCH model with Gaussian shocks and Gaussian id- 
iosyncrasies produces a slightly steeper 5-year correlation skew than the double-t 
copula, even when the latter is taken with fat-tailed idiosyncrasies. When we 
turn on the Student-t return residuals for the market factor dynamics (see figure 
[8]) the differences in the 5-year skew become quite significant. 

The explanation of the correlation skew in the dynamic TARCH-based mod- 
els is similar to the static double-t copula when one considers a particular time 
horizon. The separation of aggregate returns for the common market factor 
and idiosyncratic factors remains valid for all time horizons. The diminishing 
importance of the idiosyncratic returns compared to the market factor for the 
greater values of detachment level K explains most of the steepness of the corre- 
lation skew. Practitioners using the static models often have to assume heuristic 
term structure dependence for base correlations, typically without fundamen- 
tal reasons why one choice or another is preferred, and consequently leading to 
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biased relative value assessments between tranches of different maturities. Our 
TARCH-bascd model, in contrast, produces a characteristic pattern of correla- 
tion skew dependence on T which is driven by the speed of convergence of the 
idiosyncratic factors to normal and can serve as a starting point for such relative 
value assessments. We confine our discussion of this subject to highlighting of 
this possibility, since its detailed analysis in the context of actual market prices 
is well beyond the scope of our paper. 

4.3 Sensitivity to model parameters and hedging applica- 
tions 

The predictable shape of the correlation surface, demonstrated above, is comple- 
mented by an equally important feature of the dynamic loss generating models - 
they lead to a well-defined dependence of the correlation skew on model param- 
eters. The ability to calculate the sensitivity of expected losses to underlying 
parameters is crucial in risk management applications. In particular, sensitivity 
of the correlation surface with respect to the underlying portfolio hazard rate 
leads to a non-trivial adjustment of the CDO tranche deltas. 

First, let us demonstrate this sensitivity for TARCH model with Gaussian 
idiosyncrasies. Figure |oj shows the correlation skew p (^K ', p t (h) , Rj of the 5- 

year tranches with various detachment levels if as a function of varying portfolio 
hazard rate h. The range of variation is chosen from lOObp to 500bp, which 
corresponds roughly to portfolio spreads ranging from 40bp to 200bp, with 
R = 0.4. From the visual comparison of figures [9] and [7] it appears that the 
dependence of the correlation skew for a fixed term to maturity but varying level 
of hazard rates is very similar to the dependence of the correlation surface on 
the term to maturity. This similarity is natural, since the first order effect is the 
dependence on the level of the cumulative default probability p t (h) = 1 — e~ ht 
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Figure 8: Correlation surface for TARCH model (a=0.004, a d =0.094, /3=0.927) with 
Student-t shocks (v=8. 3) and the slices of the surface for 1, 3, 5 and 7 year maturities. 



which depends on the product of h ■ t rather than on the hazard rate or the 
term to maturity separately. For each level of this product, we get a specific 
level of the default threshold in the structural credit risk model. The higher this 
threshold, the closer is the sampled region to the center of the latent variables 
distribution and the less it is affected by the tail risk - thus leading to a lower 
level and natter skew of the credit correlation. 

However, there is a second order effect which makes these two dependencies 
somewhat different. Let us recall first that idiosyncrasies with less fat tails 
(higher degrees of freedom) correspond to flatter correlation skew, as we argued 
in the previous section. Since the idiosyncrasies converge to normal distribution 
faster than the market factor as the return aggregation horizon grows, we can 
deduce that the dependence on the term to maturity with fixed hazard rate 
should exhibit a faster flattening of the correlation surface than the dependence 
on the hazard rate with fixed term to maturity. 

The right hand side figure in [9] shows a comparison of the change in corre- 
lation when going from 5-year horizon to 10-year horizon with constant hazard 
rate set at lOObp, against the change in correlation of fixed 5-year slice when 
hazard rate goes from lOObp to 200bp. One can observe that the term-to- 
maturity extension indeed causes a greater degree of flattening than the hazard 
rate increase. 

The dependence of the correlation surface on hazard rate has a strong ef- 
fect on CDO tranche deltas. The precise calculation of tranche deltas in our 
framework requires one to derive the tranche loss probabilities from the shape 
of the correlation surface, as outlined in section [2~2} and then use the standard 
pricing techniques described in Schonbucher |27j . However, the magnitude of 
the adjustment can be estimated more easily by neglecting the interest rates and 
looking only on the protection leg of the given equity tranche. The tranche delta 
is dominated by the hazard rate sensitivity of the present value of its protection 
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Figure 9: Left figure: the dependence of the 5-year correlation skew on the level of the 
hazard rates. Right figure shows 3 correlation surface slices, contrasting the flattening 
of the skew with growing hazard rates and term to maturity. 



leg, which in turn is proportional to the expected tranche loss: 



A(o„ff] (pt) oc 



dEL 



(0,K] 



(*) 



dh 



E%L (0iK] (t)+p h [K,p t ,R\ E^L [0 , K] (t) (42) 



where, in accordance with the definition of the correlation surface, we have: 

EL {0>K] (t) = E G L [Q . K] {p (K, Pt (h) , i?) , p t (ft)) (43) 

Let us define the delta adjustment factor as the percentage adjustment to 
the tranche delta compared to fixed-correlation Gaussian delta of the tranche: 



\o,k] (Pt) = Afo, K ] (Pt) (1 + S adj (K, Pt )) 



(44) 



Since the first term in eq. (42 1 corresponds to the Gaussian delta of the 



tranche, we can see that the delta adjustment factor is equal to the correlation 
surface sensitivity times the tranche loss sensitivity ratio: 



S a dj (K,p t ) 
where 



Ph 



(k,p u r) 



E°L (0tK] (t) 



Pi, 



(K,pt,R)*{K,pt) (45) 
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E^L^ K] (t) = - l -^^^- l {p t ),-d l -~^) (46) 

3? W) - W> = (1 - P.) * (l - A) * ( - dl+ ^" 1(Pt) 

(47) 

* - -M" 1 (P.) ' ( — ) (48) 

and the functions are evaluated at p = p (K,Pt (h) , Rj ■ Thus, the tranche loss 
sensitivity ratio is given by: 



*{Kp t ) = E ° L ^ t) = 1 H^M'-^-VP) m 

* {K > Pt) E%L (0tKi (t) 2(l-p t )t^p [ ' 

Figure[lOl illustrates these calculations for the case of TARCH-based dynamic 
loss generating model. We can see that the systematic delta estimates in our 
model are anywhere from 40% to 100% greater than the conventional fixed- 
correlation delta estimates stemming from the Gaussian copula model. The 
adjustment is greatest for the lowest values of K, and drops quickly as the 
detachment level grows. While comparison with the actual tranche market price 
sensitivity is beyond the scope of our paper, we would like to note that both the 
the significant under-estimation of the equity tranche deltas and the relatively 
smaller amount of the delta error for more senior tranches are in line with the 
market experience during the correlation dislocation in May-September of 2005. 



5 Summary and Conclusions 

In this paper we have introduced and studied a new class of credit correlation 
models defined as an extension of the structural credit model where the latent 
variables follow a factor-ARCH process with asymmetric volatility dynamics. 
To build the foundation for our model, we have studied the time aggregation 
properties of the multivariate dynamic models of equity returns. We showed that 
the dynamics of equity return volatilities and correlations leads to significant 
departures from the Gaussian distribution even for horizons measured in several 
years. The asymmetry appears to "survive aggregation" longer than fat tails 
based on the parameters estimated from the real data. The main source of 
skewness and kurtosis of the return distribution for long horizons is the dynamic 
asymmetry of volatility response to return shocks. 

We introduced the notion of the correlation surface as a tool for comparing 
loss generating models, whether defined via a single-period (static) copula, or 
via multi-period (dynamic) latent-variable framework, and for simple and con- 
sistent approach to non-parametric pricing of CDO tranches. We showed that 
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Figure 10: Left figure: the dependence of the tranche loss sensitivity ratio and 
on K . Right figure shows the product of the delta adjustment factor as a function 
of K. The base case corresponds to TARCH model (a=0.004, a d =0.094, /3=0.927) 
with Gaussian idiosyncracies, aggregated over 5 year horizon. Portfolio hazard rate is 
equal to 100bps. 



portfolio loss distributions with smooth pdf can be easily reconstructed from 
the correlation surface using its level and slope along the K-dimention. 

We considered the differences in the correlation surfaces generated by static 
models, including Gaussian, Student-t and Double-t copula, and dynamic mod- 
els including GARCH(1,1) with Gaussian and Student-t shocks, and TARCH(1,1) 
with Gaussian and Student-t shocks. From this comparison we can conclude 
that the most relevant stylized facts for explanation of the market observable 
correlation skew are, in order of their importance: 

• the independence of the market factor and idiosyncratic returns (no com- 
mon mixing variables); 

• the persistent asymmetry of aggregate return distribution of the market 
factor, which in case of TARCH models occurs as a consequence of the 
asymmetric volatility dynamics; 

• the fat tails in the market factor returns; 

• the fat tails in the idiosyncratic returns; 

• the slower convergence to Normal distribution of the market factor com- 
pared to idiosyncratic returns; 

Importantly, in our dynamic framework, the correlation surface is not only 
explained, but predicted - based on empirical parameters of the TARCH process 
and the parameters describing the reference credit portfolio. The model also 
predicts a specific sensitivity of the correlation surface to changes in various 
parameters, including the average hazard rate of the underlying portfolio. 
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The inability of static models to incorporate changing base correlations are 
at the heart of the difficulties faced by these models during the credit market 
dislocations. In particular, our model reveals that the systematic deltas of equity 
tranches are understated by the industry standard static copula models, since 
the growing portfolio spread (hazard rate) should lead to an additional drop 
in equity tranche prices due to decreasing implied correlation level. Similarly, 
the static models require making additional assumptions about the correlation 
skew at different maturities. The early market convention of keeping this skew 
constant which some practitioners still adhere to is very far from realistic as can 
be seen from the results of section |4~2) The more reasonable assumption is that 
the correlation level decreases and the skew flattens for longer horizons. 

We should note however, that these conclusions are based on an implicit 
assumption that the model includes a single common return factor, and that 
the parameters of the dynamic model are constant over time. While this is a 
weaker assumption that an outright imposition of the constant correlation skew, 
it may still be too strict in some circumstances. A possible direction for gener- 
alization of our model is to move from a single market factor to a multi-factor 
framework which can make the model much more flexible in terms of both the 
detachment level K and term structure T dependence of the correlation surface 
p = p(K,T). The well-documented importance of both macro and industry 
factors for explanation of equity returns suggests that such a generalization is 
not only desirable from calibration point of view but also warranted empirically. 
While the analytical tractability of the model will suffer, the numerical accuracy 
will likely remain intact when using Monte Carlo simulations. 

Whether in a single factor or a multi-factor setting, many of our conclusions 
reflect the limitations of the large homogeneous portfolio approximation which 
we have adopted in this paper. In particular, it is clear that even determin- 
istic but heterogenous idiosyncrasies, market factor loadings and hazard rates 
could lead to significant changes in portfolio loss distribution and consequently 
to the correlation surface of the model. An extension of our model to such het- 
erogenous case is possible, although the computational efforts will increase very 
significantly. 

In conclusion we note that the ARCH family of time series models [H] had 
proven quite successful in explaining the behavior of implied volatility smile and 
skew and stock index option pricing |13j . Given the similar empirical motiva- 
tion of our model, and multi-faceted analogies with equity derivatives pricing 
outlined throughout the paper, we believe that our approach can lead to similar 
advances in the portfolio credit risk modeling, and shed new light on pricing of 
CDO tranches. 
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A Kurtosis and Skewness of Aggregated TARCH 
Returns 

T 

In this notes we analyze kurtosis and skewness of aggregated returns Rt — ^""Vt 

t=i 

when r t is assumed to follow TARCH(1,1) process 
r t = cr t s t 

of = (1 - C) v 2 + cwf-i + adri_ 1 l {rt _ 1 <o} + /3er t 2 _i 

where returns innovations St are assumed to be iid, have zero mean and unit vari- 
ance. We are interested in variance, skewness and kurtosis of time aggregated 
returns. To make sure that those moments are finite we need corresponding 
moments of the return innovations to be finite. Particularly, we assume that e t 
has finite kurtosis. Let us introduce the following notations for the central and 
truncated moments of e t 





= E(e t ) = Q 




= E(e*)=l 


< 


= E { e t l {e t <0}) 


Se 




si 


= E ( £ t 1 {et<0}) 








= E(etl {et < 0} ) 



Lemma 7 The following recursions hold for TARCH(1,1) model 
cov t _ x (r^r 2 t+u ) = pcov t _ x (rf , r 2 t+u _ x ) for u>l 
covt-i ( r t r t+i) = ctvan-! (r t fc+2 ) + a d var t ~! (r t fc+2 l { ,, t < 0} ) 

Proof. 

cov^ (r t fe , r 2 t+u ) = cov t _ x (r k t [(1 - Q ^ + ar 2 +u _ x + a^jl,,,^^^) + /3cr 2 + , 
= + acovt^ (r t fe , r t 2 +M _ 1 ) + a d cov t _ x (r t fc , r 2 +ll _ 1 l {rt+u _ 1 < 0} ) + /3ccw t _i (r t fe , of +u _ 
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if u> 1 then 



covt-x {r k t ,rl +u _ x \ rt+u _ xm ) = v d cov t -i (r t fe , r 2 t+u ^) 
covt-i (rt,crt + u-i) = covt-i (r t fe ,r t 2 +u _ 1 ) 



If u=l then 



cov t -\ (r^a^.i) = 



Proposition 8 Suppose < £ < 1 and the return innovations have finite skew- 
ness, s e , and finite "truncated" third moment, s d , then conditional third mo- 
ment of T-period aggregate return Rt,t+T has the following representation for 
TARCH(1,1) 

EtRl t+ T = s e £> (°*+«) + 3 (as £ + a D sf) ~f ] U E t (a] +u ) 

u=l u=l ^ 



In addition if Ea\ is finite then unconditional skewness of Rt t t+T is given by 



ER 



t,t+T 



E{Rlt + T? 12 



1 o 1 i d\ T (! - - 1 + P T 



Jl/2 e 1 J"3/2 



(W) 2 



E 



Proof. Using Lemma 7 we have 



t+T 



M E 



Vu=t+1 



E i I r ti r t 2 r t3 

K t+l<t 1 <t 2 <t 3 <t+T 

X> 3 + „ + £ 3 ^K<) 



u=l 
T 



t+l<ti<t 2 <t+T 



E^ + 3 E f 2 ^ 1 { aE * K) + ^ (^{ rtl <o} 

u=l t+l<ti<t 2 <t+T 

5> t (r t 3 +u ) + 3]T ^ c (r t 3 +u ) + a d E t (r 3 +M l {rt+ „< 0} )) 



u=l 



i - c T - u 



= s e J2E t (a t 3 +M ) + (as £ + a d s d £ )J2 ^T7 Et ^+u) 



u=l 



Using the law of iterated expectations 



t+T 



t+T 



El ]T r J = £ U E r - 



Vu=t+1 



Ts e + 3 (as £ + ays" 



r(i - Q - i + 
(i-C) 2 



St is iften computed using the simple formula for the unconditional variance 
E(Rl t+T )=a\ m 
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To derive unconditional kurtosis we define the following unconditional auto- 
correlations 



In = 1-n = corr(r 2 _ n ,r 2 t ) 
ip n = corr(r t - n ,r 2 ) for n > 1 
i> itj ee E (r t -ir t -jr1) for 1 < j < i 
Lemma 9 j n , ip n and ipij decay exponentially as n and i — j increase 
In = Cln-i = C" _1 7i for n > 1 
(fin = Cfn-1 = C"~Vi f or n > 1 
4>i,j = CV'i-ij-i = C J_ V»-j+i,i for 1 < j < i 
where 71, <pi and ipk,i are given by 

71 = a (k r -l)+a d (kf. - vf) + /3k r /k £ 
ipi = as r + ctdS r 
t/j k ,i = olE {r t - k +\rf} + a d E (r t _ fe+ irfl {rt < 0} ) 

with v d _ ^Mhr^y) s _ E ( r ") s d _ E ( r ^T<Q}) k _ and 

Willi U £ — Er 2 , b r — (Er 2 ) 3/2 ' r ~ (Er 2 ) 3/2 ' r ~ (Er 2 ) 2 

~ (Er 2 ) 2 ■ 

Proposition 10 If 

C = E(p + ae 2 t + a D e 2 l {£t < 0} ) = f3 + a + a D v d < 1 

f ee E (/3 + ae 2 + a D e 2 l {et < 0} ) 2 = f3 2 + a 2 k e + a 2 D k d E + 2af3 + 2a D /3v d + 2aa D k d e < 1 
then unconditional kurtosis of r t , K\, is finite and 

Proof. If the 4th moment of r t exists then the following equation must hold 
Ert = E(st)E(at) 

= k £ E ((1 - C) o- 2 + ar 2 _ 1 + adrf-^ir^Ko} + /3cr 2 -i) 2 
= k £ ((l - C) 2 o- A + 2 (1 - C) a 2 E (ar t 2 _i + a^-il^^o} + /3a 2 -i) 
+ ( ar t-i + a dr 4 2 _ 1 l {rt _ 1 < } + fio-t-i) ) 
= k £ ((1 - C) 2 a 4 + 2 (1 - C) P <J 4 + ZEaU) 
Therefore Erf necessarily solves 

Ert = k £ {l-C 2 )o- 4 + ZErt 



33 



Proposition 11 If the distribution of e t is symmetric and ad = then uncon- 
ditional kurtosis of Rt, if exists, is given by the following formula: 

Kt = 3 + - 3) + 6^ r(1 a ° ^^ for T > 1 (50) 

1 -C 2 

K 1= k £T ^ (51) 
w/iere k e is unconditional kurtosis of e t and 

£ = E(f3 + a£ 2 t + a d £ t 2 l {£t < 0} ) 2 = /3 2 + a 2 k e + a 2 d k d e + 2af3 + 2ay/3^ + 2aa d k d e . 
7l ee corr (r t 2 _ l5 r 2 ) =a(k r -l) + a d (kf - vf) + f3k r /k £ 

Proof. 

=E^(0+e E 

\u=t+l / u=l t+l<ti<t 2 <t+T 

= X> (O + e E [«* K ■ <) + * K) E «)] 

«=1 t+l<ti<t 2 <t+T 

= TE (rt) + 6^^E (r 2 ) 2 + 6cov (r^r 2 ) £ C* 2 "* 1 " 1 

t+l<ti<t 2 <t+T 

^ Ti. W ) + 6^1 E (rtf + 6cov (r^r 2 ) 
substituting the derived 4-th moment into the definition of the kurtosis Kt = 

4 



/ t+T \ * 

S ( r u J /.E (r 2 ) completes the proof. 

\u=t+l ) 



B Estimation Results for SP500 

In this appendix consider the estimation results of several TARCH(1,1) speci- 
fications for SP500 weekly returns, to provide empirical context for the rest of 
the paper. We obtained the daily levels of SP500 from CRSP database. The 
total number of observations is 10,699 and covers the period from 07/02/1962 
till 12/31/2004. We constructed weekly log returns and estimated the parame- 
ters of TARCH and GARCH models with Gaussian and Student-t shocks for 2 
samples - full and post-1990. 

Tables 1 and 2 shows estimated parameters and various data statistics. Note 
that the Student-t distribution has an additional parameter, degrees of freedom 
v, that adjusts the tails of the error distribution. Since the Gaussian distribution 
is nested within the Studcnt-t as a limit of large degrees of freedom, and since 
the estimates of the full unconstrained model result in a relatively small and 
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statistically significant value of the degrees of freedom, we conclude that the 
data points toward the fat-tailed return shock distribution. 

On the other hand, the asymmetric TARCH model is nested within the 
symmetric GARCH in the limiting case ad — 0. The estimated asymmetric 
coefficient ad in the TARCH model is not only non-zero, but significantly higher 
than the symmetric coefficient a for both complete and post 1990 samples, both 
daily and weekly frequencies and Gaussian and Student-t shock distributions. 
Thus, we conclude that the asymmetric volatility is prominently present in the 
data. The best fit model among those considered is the TARCH(1,1) with 
Student-t distribution of return innovations. The additional parameters of this 
model are statistically significant. 

To make sure that asymmetry in volatility is not a result of several extreme 
negative returns like 1987 crash we provide data statistics and re-estimated 
parameters of TARCH models for trimmed full and post 1990 samples. The 
trimming is done by cutting excess volatility in the most extreme 0.1% obser- 
vations of both positive and negative return. 



Table 1: SP500 moments. 



Sample period 


Daily 


Weekly 


1962-2004 
1990-2004 


Sf S r ky Vy. 
-1.40 -2.43 39.83 0.53 
-0.11 -1.14 6.67 0.51 


Sy Sy, ky xi y, 
-0.55 -1.35 7.01 0.55 
-0.64 -1.36 6.10 0.56 


SP500 moments( After trimming 0.1% of extreme positive and negative returns ) 


1962-2004 
1990-2004 


0.05 -1.03 5.95 0.50 
0.04 -1.01 5.56 0.50 


-0.39 -1.18 5.26 0.54 
-0.50 -1.24 5.22 0.55 



Table 2. 



Estimated parameters of GARCH(1,1)/TARCH(1,1) with Gaussian/Student-t shocks 
on weekly SP500 returns. The total number of return observations is 756 for post-1990 
sample and 2,139 for the full sample starting in 1962.. Number below each parameter 
estimate in parenthesis is an asymptotic standard deviation, LogL is corresponding 
loglikelihood value. GARCH parameters correspond to the volatility specification: 
o\ = oj + ar\_ x + adrt_ l l{ rt l <Q} + j3a^_ l . vis degrees of freedom of Student-t 
distributed return innovation e t . The constant term of the volatility process is not 
shown since it is not used in the simulations. 
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Sample Dates 


01/01/1990-12/31/2004 


01/01/1962-12/31/2004 


Model 


GARCH 


TARCH 


TARCH + t 


GARCH 


TARCH 


TARCH + t 


a 


0.044 


0.007 


0.004 


0.107 


0.037 


0.032 




(0.0073) 


(0.024) 


(0.019) 


(0.013) 


(0.016) 


(0.0125) 


Old 




0.112 


0.094 




0.136 


0.106 






(0.046) 


(0.033) 




(0.033) 


(0.0223) 




0.953 


0.918 


0.927 


0.886 


0.877 


0.894 




(0.0003) 


(0.0022) 


(0.0013) 


(0.0022) 


(0.0031) 


(0.0015) 


V 






8.31 






10.19 








(2.75) 






(2.42) 


LogL 


1855.6 


1861.7 


1877.4 


5347.4 


5368.6 


5397.7 



C Monte Carlo Simulations 

Most of the numerical estimates for credit risk in this paper are obtained by 
Monte Carlo simulation. Here we outline the simulation procedure for two such 
calculations, the estimation of the pairwise default correlation coefficient, and 
the estimation of the tranche losses under the LHP assumption. 

The default correlation coefficient, p d (p) for the factor GARCH and TARCH 
models is calculated based on the simulated factor time series and closed form 
formulas of conditional default probabilities: 

• simulate the common factor, R m T , I = 10,000 times and normalize it to 



have variance 1 



for each p find d T (p) that solves 7 ( ^ T /rzri' T ) = P 



=i 



• calculate p d (p) = where p 12 = ^ dr(p) bR ™- T 

i=l 

The likelihood bounds in Figure[3]were obtained by repeating this procedure 
1000 times. 

When simulating portfolio loss distribution under the large homogeneous 
portfolio (LHP) assumption we again begin by simulating the aggregated market 
factor return. The latent variables are assumed to have symmetric one factor 
structure with the factor following TARCH (1,1) model. For each realization of 



the market factor the portfolio loss is given by the LHP formula ( 10 ) 

(!-*)*( 



Lt = 

where 



T / T 



Rm.T — ^^frn,u/ std I yY m ,s ) is a normalised return over horizon T 

«=i \s=l / 

generated using time aggregation of simulated TARCH (1,1) returns with 

unconditional volatility equal to 1 
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• d,T is calibrated so that the probability of J?j ; T = 6i? m ,T + Vl — & 2 £r 
hitting is equal to single name default probability 



• b is the factor loading that is chosen to match a given unconditional linear 
correlation p = b 2 

To calculate the expected tranche losses generated by the model and to 
calibrate d T we use I = 100,000 independent Monte Carlo simulations of the 
factor and then use corresponding sample moments: 






7E/(o,*]((i-*) 



( 




) 



) 



Vl-b 2 
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